Self-accelerating the normal DGP branch 
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We propose a generalised induced gravity brane-world model where the brane action contains an 
arbitrary f(R) term, R being the scalar curvature of the brane. We show that the effect of the f(R) 
term on the dynamics of a homogeneous and isotropic brane is twofold: (i) an evolving induced 
gravity parameter and (ii) a shift on the energy density of the brane. This new shift term, which 
is absent on the Dvali, Gabadadze and Porrati (DGP) model, plays a crucial role to self-accelerate 
the generalised normal DGP branch of our model. We analyse as well the stability of de Sitter self- 
accelerating solutions under homogeneous perturbations and compare our results with the standard 
4-dimensional one. Finally, we obtain power law solutions which either correspond to conventional 
acceleration or super-acceleration of the brane. In the latter case, no phantom matter is invoked on 
the brane nor in the bulk. 

PACS numbers: 



I. INTRODUCTION 

Understanding the recent acceleration of the universe 
is a challenging task facing the cosmological and particle 
physics community. The first evidence for the accelera- 
tion of the universe was provided by the analysis of the 
Hubble diagram of SNe la a decade ago [l[ . This discov- 
ery, together with (i) the measurement of the fluctuations 
in the cosmic microwave background radiation (CMB) 
which implied that the universe is (quasi) spatially flat 
and (ii) that the amount of matter which clusters grav- 
itationally is much less than the critical energy density, 
implied the existence of a dark energy component that 
drives the late-time acceleration of the universe. Subse- 
quent precision measurements of the CMB anisotropy by 
WMAP and the power spectrum ofgalaxy clustering 
by the 2dFGRS and SDSS surveys 0, gf have confirmed 
this discovery. 

A possible approach to describing the late-time accel- 
eration of the universe is to consider a modification of 
gravity, such that a weakening of this interaction on the 
appropriate scales induces the recent speed up of the uni- 
verse (cf. Refs. 0,0,0, [I])- In other words, the weakening 
of gravity on large scales would provide an effective neg- 
ative pressure that would fuel the late-time acceleration 
of the universe. 

The simplest of these models is the Dvali, Gabadadze 
and Porrati (DGP) scenario 0, [l(| [Hi, which corre- 
sponds to a five-dimensional (5D) model. In this model, 
our universe is a brane; i.e. a 4D hyper-surface, em- 
bedded in a Minkowski space-time. Matter is trapped 
on the brane and only gravity experiences the full bulk. 
The DGP model has two branches of solutions: the 
self- accelerating branch and the normal one. The self- 
accelerating brane, as its name suggests, speeds up at 



late-time without invoking any unknown dark energy 
component. On the other hand, the normal branch re- 
quires a dark energy component to accommodate the cur- 
rent observations [u3.[l3| . Despite the nice features of the 
self- accelerating DGP branch, it suffers from serious the- 
oretical problems like the ghost issue [l4[ . The main aim 
of this paper is to consider a mechanism to self-accelerate 
the normal branch which is known to be free from the 
ghost issue 
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This mechanism will be based on a modified Hilbcrt- 
Einstcin action on the brane and the simplest gravita- 
tional option is an f(R) term. Extended theories of grav- 
ity based on 4D f(R) scenarios have gathered a lot of 
attention in the last years (cf. the extensive lists of refer- 
ences in 0,0, 0]). It has been shown that these 4D models 
should follow closely the expansion of a LCDM universe 
[l5l . lit! Qj} and could have distinctive signatures on the 
large scale structure of the universe [HI, [l!| . On the other 
hand, several methods have been invoked to reconstruct 
the shape of f(R) from observations [13, HH, HH , for ex- 
ample, by using the dependence of the Hubble parameter 
with redshift which can be retrieved from astrophysical 
observations. We will show that an f(R) term on the 
brane action can induce naturally self-acceleration on the 
normal DGP branch. 

This paper is outlined as follows. In section HU we 
present our model and deduce the modified Einstein 
equation on the brane. In section IIII1 we describe the 
dynamic of a homogeneous and isotropic brane on the 
framework introduced previously. Then we describe the 
effect of an f(R) term on a standard induced gravity 
brane. On the next section lTVl we obtain self-accelerating 
solutions corresponding to de Sitter space-times for both 
generalised DGP branches; i.e. the generalised normal 
solution (now sclf-accclcrating) and the generalised self- 
accelerating solution. We study the stability of these so- 
lutions under homogeneous perturbations. We compare 
also our results with the standard 4D ones. In section 
IVl we construct power law solutions for the brane ex- 
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pansion. Finally, in the last section we summarise and 
conclude. 



Here, Ti^ v corresponds to the quadratic energy momen- 
tum tensor [241 



II. A GENERALISED INDUCED GRAVITY 
SCENARIO 
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We consider a brane, described by a 4D hyper-surface 
(h, metric g), embedded in a 5D bulk space-time (B, 
metric g^), whose action is given by 



S = 



where k| is the 5D gravitational constant, i?[g( 5 )] is the 
scalar curvature in the bulk and K the extrinsic curvature 
of the brane in the higher dimensional bulk, correspond- 
ing to the York-Gibbons-Hawking boundary term [23|. 
For simplicity, we will assume that the bulk contains 
only a cosmological constant; i.e. £5 = —U. Therefore, 
the bulk space-time geometry is described by an Einstein 
space-time 
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The 4D Lagrangian corresponds to 
£4 = ctf(R) + C m , 



(2.2) 
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where R is the scalar curvature of the induced metric 
on the brane, g, and a is a constant that measures the 
strength of the generalised induced gravity term f{R) 
and has mass square units. Notice that therefore the 
function f(R) has mass square units. On the other hand, 
C m corresponds to the matter Lagrangian of the brane 
which in particular may include a brane tension. The pre- 
vious action, includes as a particular case the DGP model 
0, [3 when the bulk is flat, f(R) = R and a = \/2k\ 
where k| is proportional to the 4D gravitational con- 
stant. For the time being we will consider / an arbitrary 
function of the scalar curvature of the brane. 

As we are mainly interested in the cosmology of a ho- 
mogeneous and isotropic brane it is quite useful to fol- 
low the approach introduced by Shiromizu, Maeda and 
Sasaki in 1 [24|. Then, the projected Einstein equation on 
the brane reads, where we have assumed a Z^-symmctry 
across the brane, 



and E^ v is the (trace-free) projected Weyl tensor on the 
brane. 

The total energy momentum on the brane is defined as 
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and can be split into two terms 
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The first term Tj^ 1 corresponds to the energy momentum 
tensor of matter (which include in particular the brane 
tension) on the brane. The second term 
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corresponds to the energy momentum tensor due to the 
generalised induced gravity term, f(R), on the brane. 

Now, if / is proportional to the scalar curvature of the 

( f ) 

brane, then t^J is proportional to the Einstein tensor of 
the brane; i.e. the standard induced gravity brane-world 
scenario is retrieved: 



(2.9) 



Using the 5D Codacci equation, the bulk Einstein 
equation, and the junction condition at the brane, it 
turns out that the total energy momentum tensor of the 
brane is conserved t^ v (24|, i.e. 



vv = 0. 

On the other hand, because 2 

VtJJ? = 0, 



(2.10) 



(2.11) 



we can then conclude that the energy momentum tensor 
of matter on the brane is conserved 



V^ T ( m ) = 



(2.12) 



1 



(2.4) 



1 For an alternative approach to deduce the equations of evolution 
of a DGP brane with curvature modifications on the brane action 
see [25ll2fj|. See also [2?} for a brane-world model with an f(R) 
term. 



2 We have proven this equation by making use of the 4D Bianchi 
identity on the brane; i.e. VG M „ = 0, and the relation between 
the non commutative character of two covariant derivatives and 
its relation to the Riemann curvature tensor (again on the brane), 
see for example equation 3.2.12 of Ref. [28ll , Therefore, the con- 
servation relation (|2.11| l can be proven in analogy to how it is 
done in the standard 4D f(R) scenario. 
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III. DYNAMICS OF A FLRW BRANE 

In what follows, we consider a Friedmann-Lemaitre- 
Robcrtson- Walker (FLRW) brane. If the brane is ho- 
mogeneous and isotropic it is known that the most gen- 
eral vacuum bulk in which the brane can be embedded 
is a Schwarzschild anti-de Sitter space-time p9j . In this 
particular case, it is also known that the tensor E^ v is 
conserved and by virtue of its traceless property, the pro- 
jected Weyl tensor on the brane behaves as a "dark" ra- 
diation fluid. 

The matter sector on the brane, prescribed by , 
can be described by a perfect fluid with energy density 
p( m ) an d pressure p( m > , where p^ is conserved as a con- 
sequence of Eq. (|2.1ip . On the other hand, an effective 

energy density and an effective pressure associated to t^J 
can be defined as follows 
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where the label i makes reference to the spatial coordi- 
nates of the brane. The parameter k = ±1,0 depending 
on the geometry of the spatial sections of the brane. The 
dot refers to a derivative respect to the cosmic time of 
the brane t and the prime refers to a derivative respect 
to the scalar curvature R. 

We recover the evolution equations of a homogeneous 
and isotropic universe in the framework of 4D f(R) theo- 
ries (given for example in the recent review Q) by simply 
imposing 
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This is a simple consequence of the definition (|2.6p . In 
the standard 4D case the tensor (|2.6[) has to vanish due 
to the principle of least action because £4 would be the 
total Lagrangian and therefore the variation of the 4D 
action with respect to the metric would have to be zero. 

In the brane-world scenario things get a bit more in- 
volved because of the presence of (i) the quadratic energy 
momentum tensor LI^ and (ii) the projected Weyl ten- 
sor E^ v in the modified Einstein equation on the brane 
(|2.4[) . As we have already mentioned the effect of the 
projected Weyl tensor for a FLRW brane (embedded in a 
Schwarzschild anti-de Sitter space-time) will be the pres- 
ence of a "dark" radiation fluid in the evolution equations. 
On the other hand, the quadratic energy momentum ten- 
sor LT^j, contributes to the effective Einstein equation 
(O with 
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In the previous equations 
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(3.5) 



where p^> and p^ are given in Eqs. (|3.ip and (|3.2p . 
So, finally we can write down the modified Friedmann 
equation on the brane as 



H 2 + \ 
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The second term on the right hand side (rhs) of equation 
(|3.6j) corresponds to the "dark" radiation energy density. 
The spatial component of Einstein equation can be ex- 
pressed as 



2H + 3H 2 



C 
3aA 



12 



p(p + 2p), (3.7) 



where the energy density p and the pressure p are de- 
fined in Eq. (|3~5]) . Even though Eqs. ([3~6]) and (|3~7|) look 
very simple this is not case. For example, the brane Ray- 
chaudhuri equation (|3.7p contains the effective quantities 

In order to compare this model with the standard in- 
duced gravity brane-world scenario [3(| , i-e. f(R) = R 
in the Lagrangian (|2.3p . it is useful to define a rescaled 
induced gravity coupling 

a (r) = af (3.8) 

and split further the effective quantities p^> and p^ as 



where 
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The magnitudes p^ and p^ , we just defined, get re- 
duced to the effective energy density and pressure associ- 
ated to the standard induced gravity models; i.e. models 
with f(R) = R. The effect of an f(R) term on an in- 
duced gravity brane is twofold: (i) an evolving induced 
gravity parameter, a^ r \ and (ii) a shift on the matter en- 
ergy density and matter pressure quantified by p^ and 
p( c \ respectively. In fact, we can rewrite the Friedmann 
equation as 
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H 2 
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(3.14) 



The previous expression shows the existence of two branches of solution for H 2 as a function of the effective energy 
density of the brane p( m ' + p^ . On the other hand, the modified Raychaudhuri equation is 
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j(m) + (c) _ g a (r) / #2 + 
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Thus the modified Einstein equations can be written sim- 
ilarly to that of a standard induced gravity brane [30| . 
The difference of course is hidden in i.e. an evolving 
induced gravity parameter, p^ and p^ c \ which define an 
effective energy density and pressure on the brane: 

p eS = p (ra) + p (c) ^ pC ff = p (m) + p (c) _ (3.16) 

A similar situation happens in an induced gravity sce- 
nario with an a non-minimally coupled scalar field on 
the brane [3lj . 

A very important issue that we have not yet discussed 
is the relation between the effective 4D gravitational cou- 
pling and the 5D gravitational constant k 2 . This issue 
depends strongly on which regime we are considering on 
the brane: early universe or late-universe (cf. for example 
(IcL l3lL I32I] ) . The approach used in [3l[ can be extended 
to our model to retrieve the 4D effective gravitational 
constant at early-time on the brane. Please, notice in 
that case one is considering a non vanishing brane ten- 
sion. Here we are rather interested in the late-time evolu- 
tion of the brane. Then, in this case, the brane effective 
gravitational constant can be easily obtained following 
the standard approach in the DGP model [l(|. For sim- 
plicity, we choose k = 0, U = and C = 0, then the 
modified Fricdmann equation (|3.6[) can be rewritten as 

H 2 = f p^ + pV) ± -^H. (3.17) 

Therefore, the effective gravitational constant on the 
brane is 87rG c ff = l/{2a^) which is rescaled by an 
f'(R) term with respect to the standard DGP model. 
On the other hand, we can define as well a new gener- 

(r) 

alised crossover scale, r y c , which splits the 5D regime 
from the 4D regime 3 , such that = K^a^h There- 

(r) 

fore, it is this quantity, r c , that relates the gravitational 
constant on the brane with the fundamental 5D quantity 



3 Here by a 5D regime we are referring to H oc p cff ; i.e. r£ H, 
while in the 4D regime H 2 oc p cB ; i.e. H < r [ c r >. 



«§. The precise shape and magnitude of 7c has to be 

(r) 

fixed by observation. Notice that r c ; is related to the 
crossover scale, r c , of the DGP model by = f'(R)r c . 
On the other hand, the extra-dimension includes a new 
term (±l/(ri r ^)H) in the modified Friedmann equation 
compared to 4D f(R) models (see Eq. (|3.17[) ). 

Before concluding, we point out that the two branches 
of the DGP model H,[l(j are particular cases of Eq. (|3.14p 
(and of course also of Eq. (|3.17p ). These solutions corre- 
spond to a = l/(2/c|), f(R) = R; therefore p^ and 
would vanish, C = and U = 0. These solutions in the 
absence of matter correspond to a de Sitter solution (self- 
accelerating branch) or a Minkowski space-time (normal 
branch). We will next show that the normal branch can 
become self- accelerating due to the f(R) contribution on 
the brane action. 



IV. DE SITTER BRANES 

One of the most puzzling problems nowadays in physics 
is the issue of the late-time acceleration of the universe. 
A possible approach to tackle this problem is within the 
frame-work of self-accelerating universes; i.e. could it be 
that a modification of gravity at late-time and on large 
scale be the cause of the current inflationary phase of the 
universe? A de Sitter universe is the simplest cosmologi- 
cal solution that exhibits acceleration and therefore it is 
worthwhile to prove the existence of this solution in our 
model and study its stability. This would be a first step 
towards describing in a realistic way the late-time acceler- 
ation of the universe in an f(R) brane-world model. This 
approach will also enable us to look for self-accelerating 
solutions on the modified normal DGP branch (— sign in 
Eq. (EH). 

In this section, we first obtain the fixed points of the 
model corresponding to a de Sitter space-time and then 
we study their stability under homogeneous perturba- 
tions. For simplicity, we will consider the spatially flat 
chart of the brane; i.e. k = 0, and no dark radiation 
on the brane; i.e. the bulk corresponds to a 5D maxi- 
mally symmetric space-time. Notice that even in more 
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general cases the dark radiation term will have no influ- 
ence on the late-time dynamics of the brane as this term 
is constrained to be already subdominant by the time of 
nucleosynthesis [H|. 



A. Background solutions 

The Fricdmann equation (|3.6p implies that if the brane 
geometry corresponds to a de Sitter space-time with Hub- 
ble rate Hq, then 
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2a Q/ - 3H 2 /'^ = Constant. 



(4.1) 



In particular, the matter energy density has to be con- 
stant. Now imposing the conservation of the matter con- 
tent of the brane, it turns out that matter in this case be- 



haves like a cosmological constant p^ r 
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3< m ). A con- 



stant matter energy density can always be re-absorbed in 
the f(R) related terms. For this reason, we will disregard 
the matter content in our analysis of de Sitter branes. Fi- 
nally, using Eq. (|3.14| we get that the Hubble parameter 
can be expressed as 



24a z FjH£ = 1 + -4a z F (R a F Q - /„) 



\4aF [a(RoF - /„) - 4aF U] 



(4.2) 



where e = ±1, the subscript stands for quantities 
evaluated at the de Sitter space-time, i?o = 12i?o an d 
F = df/dR. We recover the DGP model for f(R) = R 
and U = 0. In fact, in that case, the de Sitter self- 
accelerating DGP branch is obtained for e — 1 and the 
normal DGP branch or the non-self-accelerating solution 
for e = — 1. When the brane action contains curvature 
corrections to the Hilbcrt-Einstcin action given by the 
brane scalar curvature, the branch with e = — 1 is no 
longer flat and accelerates (cf. Figs [1] and \2§ . Therefore, 
an f(R) term on the brane action induce in a natural 
way self-acceleration on the normal branch. Most impor- 
tantly, it is known that such a branch is free from the 
ghost problem (see 14[ and references therein). The rea- 
son behind the self-acceleration of the generalised normal 
brane is twofold: (i) p(°> defined in Eq. (|3 . 1 2[) does not 
vanishes (in general) as it reduces to 
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FIG. 1: The figure shows the behaviour of the rescaled 
squared Hubble rate Ik^c^Fq Hq for the two branches that 
generalise the DGP solution versus the rescaled energy den- 
sity defined as ^nfa 2 Fq(RqFq — fa). The blue star cor- 
responds to the normal DGP branch which is flat. The red 
star corresponds to the self-accelerating DGP branch. On the 
other hand, the blue curve corresponds to the generalised self- 
accelerating branch, while the red curve corresponds to the 
generalised normal branch. 



tion it is useful to rewrite the Hubble rate (14.21) as 



Hi 



H (4) 



1-^1 + |a 2 K|F (/ - KjUF ) 



2o?4F$ 



(4.4) 



where we have substituted Rq 
is defined as 



12H 2 in Eq. and 



6F n 



(4.5) 



The solution (|4.4|) contains both DGP branches 4 ; i.e. so- 
lutions with f(R) = R, U = and H% = 0, 1/(k|« 2 ). 



The first term on the rhs of Eq. (|4.4[) , , corresponds 
to the Hubble rate of de Sitter universes in 4D modi- 
fied theories of gravity of the kind f(R) (see for example 
[13, HE])- Therefore, the presence of the extra dimension 
implies a shift on the Hubble rate (cf. Eq. (|4.4p ). 

We notice as well that the de Sitter branes are close 
to the standard 4D regime as long as H { 
implies 



iJ 2 4 -p which 



3 1- Jl + i4a 2 F (f - ^UFo) 



4 a2F ofo 



< 1. (4.6) 



and (ii) p^ c > enters the modified Friedmann equation This relation will be helpful for the next ana i ysis . 
P-14p as an effective energy density. 

Once we have seen the effect of an f(R) term on the 

DGP model, we will ask the question the other way 

around: What is the effect of an extra-dimension on a 4 Noticc as well that the existonce of two different branches is 

4D f(R) model? In order to answer the previous ques- hidden on Eq. JO). 
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FIG. 2: This a zoom of the normal branch as it appears in 
figure Q] 



B. Stability analysis 

We next analyse the stability of de Sitter solutions un- 
der homogeneous perturbations up to first order on 5 



SH = H(t) - H . 



(4.7) 



We will follow the approach used in [34[ . The perturbed 
Friedmann equation (|3.6[) implies an evolution equation 
for SH: 

SH + 'SHoSH + m 2 cS SH = 0, (4.8) 



where the square of the effective mass rn 2 s is given by 
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In the previous equation fjm = d 2 f/dR 2 , evaluated at 
the de Sitter background solution. 

We restrict our analysis to H > as for H < there 
is always an exponentially growing mode for SH implying 
that the solution H = Ho < is unstable. Then, a de 
Sitter solution with H > is stable as long as m 2 ^ is 
positive. 

It is quite useful to rewrite rn 2 ^ by substituting 
Eq. (|4~4|) into Eq. (|4~9)) . We then get 



m eS = m (4) + m i 



hift 



"pert ' 



(4.10) 



5 We do not intend to analyse the ghost issue present in the self- 
accelerating DGP model in the present paper as it is beyond the 
scope of the present paper. 



where 
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a?4F 2 
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^l + ^4F (fo-4UF ) 
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1- \/l + -a24F (f - k 2 UF ) 



The m?j^is the analogous quantity to rn 2 s in a 4D f(R) 



model (341 ] . Therefore, the extra-dimension induces a 
shift on m 2 s caused by two effects: (i) a purely back- 
ground effect due to the shift on the Hubble parameter 
(see the second term on the rhs of Eqs. (|4.4[) and (|4.9p ) 
encoded on m 2 ack and (ii) a purely perturbative extra- 
dimensional effect described by mp ert . 

In the remaining of this subsection, wc will choose 
Fq > 0, so the induced gravity parameter is positive (see 
Eq. (|3.8|) V The last supposition guarantees a positive ef- 
fective gravitational constant on the brane at late-time 
(cf. Eq. (|3.17p ). Finally, we will also assume that wc 
are slightly perturbing the Hilbcrt-Einstein action of the 
brane, i.e. /o ~ Ro- Therefore, /q is positive because 
R = 12Hq. For simplicity, we will assume as well that 
U = 0. These three suppositions (U = 0, fo, F > 0) im- 
plies that the 4D regime (cf. Eq. (|4.6[) ) is reached when 
the inequality 



1 « 4a 2 F f 



(4.12) 



is fulfilled. 



The previous inequality implies that to? k > and 



pert 



< 0. Consequently, the shift on the brane Hub- 
ble parameter respect to the standard 4D case tends to 
make the perturbation heavier; i.e. the de Sitter universe 
would be more stable. However, the pertubative effect 



encoded on 



'pert 



would make the perturbations lighter 



and therefore the de Sitter space-time would be less sta- 
ble than in the pure 4D case. By imposing the condition 
(|4~T^|) on Eqs. P~TU)) and (|4~PIj) . it can be shown that 
the extra-dimension has a benigner effect in the 4D f(R) 
model; i.e. m 2 s > m 2 ^ , as long as 

Fo < 4 /o/flii- (4-13) 

So far, we have described the effect of the extra- 
dimension on the stability of a de Sitter solution under 
homogeneous perturbation in a 4D f(R) model. It is, 
however, not possible to analyse the effect of an f(R) 
term (on the brane action) on the stability of the de Sit- 
ter self-accelerating DGP solution. This is simply due to 
the fact that the f(R) contribution carries an extra degree 
of freedom which is absent in the original DGP model or 
even GR models [HI, [lj| . That is the reason why m 2 s is 
not well defined for f(R) = R. The term vanishes 
in that case. This is not surprising as something similar 
happens in General Relativity (GR) and 4D f(R) models. 
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Indeed, m 2 A y the equivalent ra 2 ff in the 4D f(R) scenario, 
is also not well defined for f(R) = R in the standard 4D 
case which only means that the analysis docs not apply 
to GR and not that de Sitter space-time is not stable in 
GR 6 . 

Finally, let us point out that the positiveness of m 2 s 
depends on which branch we are considering, as the Hub- 
ble rate Hq depends on the specific branch [see Eq. (|4.2|) ]. 
So, a given self-accelerating solution (|4.2p with Hq > is 
stable under homogeneous perturbations as long as m 2 ff 
defined in Eq. ()4.9p is positive. 



V. POWER LAW EXPANSION ON THE BRANE 

In this section we construct power law solutions for the 
brane expansion; i.e. 



( j- 



where a ,t ,(3 ^ are constants 7 and t is the cosmic 
time of the brane. For simplicity, we will consider that 



/ = foR n , fo, n — Constants. 



(5.2) 



These ansiitze have been previously used in standard 4D 
f(R) models (H,^- 

Then the effective curvature energy density p(*> and 
pressure pv> , defined in Eqs. ()3.1|) and ()3.2|) . satisfy 



,(/) 



6[/?(n-2) + (2ra-l)(n-l)] 



a/3f 
nt 2 ' 



p if) = 2(2n-3/3)[/3(n-2) + (2n-l)(n-l)] 



(5.3) 

of 
nt 2 ' 
(5.4) 



where 



/' = nf 



6/3(2/3-1) 
t 2 



rc-1 



An equation of state parameter can be defined as 

. 2n 



w f 



(/) 1+ 3/3' 



(5.5) 



(5.6) 



This parameter is not to be confused with the effective 
equation of state parameter w c g defined through 



H 



-p a 



-3(1+U) e ff) 



(5.7) 



6 We arc grateful to V. Faraoni for pointing out this to us. 

7 We henceforth disregard the case (3 = as it corresponds to a 
static brane which is not interesting from a cosmological point 
of view. 



where po is a constant. For a power law expansion w e s 
reads 



w e s = -1 



3/3' 



(5. 



In standard 4D f(R) models the quantities p^> and 
p(f> vanish (see Eq. (|3.3[) ) impl yin g a set of constraints 



in the parameters /3 and n [36j, |37j . In the brane- world 
scenario these constraints are different as a consequence 
of the modified Einstein equations (see Eqs. (|3.6p , (|3.7p ). 
Assuming that the brane is spatially flat, i.e. k = 0, 
and the bulk corresponds to a Minkowski space-time, i.e. 
U = 0, the Fricdmann equation implies 



[6/3(2/3 - l)] 2(n_1) 



(5.9) 



(5.1) while the Raychaudhuri equation imposes the constraint 



0(3/? -2) _ 



t 2 



= «fa 2 /o/3(3/3-4n) 

x [f3(n - 2) + (2n - l)(n - 1)Y 

[6/3(2/3- i)] 2 ^'- 1 ) 

x ~. . 

t 

From the previous two equations it turns out that 
1 



/3 



16 - 3n 2 a 2 f 2 



(5.10) 

(5.11) 
(5.12) 



Then, we can conclude that 1 + w c s — 2(1 + Wf) (cf. 
Eqs. (|5.6p - (|5.8p ). The prcfactor two on the rhs of the 
previous equation is a brane effect; i.e. the energy density 
p(f> appears quadratically on the modified Friedmann 
equation. 

The brane accelerates if (i) (3 > 1 or (ii) (3 < 0. Case 
(i) corresponds to conventional acceleration and it will 
hold if 



1 < i^« 2 /o 2 < 2- 



(5.13) 



Now, if the induced gravity parameter is defined roughly 
as in the DGP model i.e. a ~ 1/(2/^) = M p, we 
obtain bounds on the allowed set of fo values for the 
brane being accelerating 



16 Ml ~ 2 32 Ml 
3 M% Jo 3 Mp ' 



(5.14) 



This bound is related to the ratio of the fundamental 
Planck mass, M5 , and the effective 4D Planck mass Mp. 
Conversely, case (ii), which is given by 



2 < i^« 2 /o, 



(5.15) 
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describes super-acceleration. By assuming that the in- 
duced gravity parameter is defined as in the DGP model 
3, we get 



32 Mf ~ 2 
3 Mp J ° ' 



(5.16) 



i.e. we obtain a lower bound for /r> 

In summary, we have shown that the brane can follow 
a power law expansion such that the brane is acceler- 
ating (1 < (3) or even super-accelerating {(3 < 0). In 
the latter case, the brane faces a big rip singularity in 
its future 38]. This happens at t = where the scale 
factor, H, H, H and H blow up. For (3 < 0, we have 
chosen the cosmic time to be negative so that the brane 
expands as time pass by. Unlike in a standard relativis- 
tic framework, in f{R) models up to the third cosmic 
time derivative of the Hubble rate appear in the evolu- 
tion equation of the universe. Consequently, it makes 
sense to speak about divergences of H and its deriva- 
tive. The super-acceleration of the brane is accompanied 
by a phantom-like behaviour; i.e. 1 + w e ff < 0. Even 
the effective fluid {p^\p^) mimics a phantom fluid; i.e. 
1 + Wf < 0, this mimicry is exclusively due to curvature 
effects and takes place although no real phantom matter 
has been considered in the model. 

Alternative ways to get a phantom-like behaviour on 
the brane are based on a screening of the cosmological 



constant 
the bulk 




13|, [39( or a flow of energy from the brane to 
In the model presented here, the phantom 
mimicry is a pure gravitational effect and it does not 
involve any of the two effects mentioned above. 



VI. CONCLUSIONS 

In this paper we present a mechanism to self-accelerate 
the normal DGP branch which unlike the original self- 
accelerating DGP branch is known to be free from the 
ghost problem 8 . The mechanism is based in including 
curvature modifications on the brane action. For sim- 
plicity, we choose those terms to correspond to an f(R) 
contribution, which in addition is known to be the only 
higher order gravity theories that avoid the so called Os- 
trogradski instability in 4D models @- 

We obtain the effective Einstein equation on the brane 
following the approach introduced by Shiromizu, Maeda 
and Sasaki in Ref. (24[. We then describe the dynamic 
of a FLRW brane and identify the two branches of so- 
lutions of the generalised induced gravity brane-world 
model. The f(R) term on the brane action results on 



Notice as well that by embedding the DGP model in a higher 
dimensional space-time, the ghost issue present in the original 
DGP model may be cured |4]J while preserving the existence of 
a self-accelerating solution |42H . See also Ref. [43l , U-ll 



an evolving effective induced gravity parameter; i.e. a 
variable gravitational constant. The same term also in- 
duces a shift on the energy density of the brane through 
the new magnitude p (c) (cf. Eqs. ([3321 and ^J^). In 
the DGP model p(°> vanishes, however, for f(R) ^ it is 
generally non zero. This is precisely the reason for get- 
ting self-acceleration in the branch that generalised the 
normal DGP solution (see Eq. (|4.2p b 

We obtain all the de Sitter self-accelerating solutions 
and study their stability under homogeneous perturba- 
tions. It turns out that the self-accelerating solutions are 
stable as long as the parameter m^ ff , defined in Eq. (|4.9p 
and related to the effective mass of the perturbation, is 
positive. This parameter can be splitted into three differ- 
ent terms: (i) the standard 4D contribution m? 4 % (see e.g. 
Ref. [3]), a contribution from a purely background ori- 
gin ?™^ ack and another one of pcrturbative origin m^ cl . t , 
all defined in Eqs. (|4.10p and (|4.1ip . For those solutions 
close to the standard 4D regime, i.e. those satisfying the 
inequality (|4.12p . TOp ert tends to make the homogeneous 
perturbation heavier (njp 0rt > 0), while rn^ adk tends to 
make the same perturbation lighter (m^ ack < 0). More- 
over, we have shown that the extra-dimension has a be- 
nigner effect in a 4D f(R) model; i.e. rn 2 s > m? 4 ), as 
long as the inequality (|4.13|) holds. 

On the other hand, we have obtained power law solu- 
tions for the brane expansion which corresponds to con- 
ventional acceleration or super-acceleration. The super- 
acceleration is achieved without invoking any phantom 
matter on the brane or the bulk. 

Last but not least, it is know that 4D f(R) models are 
not free from theoretical problems (cf. Ref. [45| for a re- 
cent account on the subject), so in constructing an f(R) 
brane-world model, we should of course try to avoid these 
theoretical troubles. We have just undertaken a first step 
towards constructing realistic self-accelerating solutions 
in the normal DGP branch. There are still many issues to 
be addressed, for example which f(R) should we pick up 
to be in agreement with the cosmological observations 
and the solar system tests? We leave these interesting 
issues for future works. 
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